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Abstract 

The radiative decay ujj — > + 7 of massive neutrinos is analyzed in the framework 
of the standard model with lepton mixing for very strong magnetic fields B 3> B CT = 
m\je ~ 4.14 x 10 13 G. The analysis is based on the approximate decay amplitude 
obtained by Gvozdev et al. Numerical results as well as analytical approximations 
for the decay rate are obtained for energies of the initial neutrino below and above 
the electron-positron pair creation threshold 2m e . 
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1 Introduction 

In the Standard Model of electroweak interactions (SM) the neutrino masses 
are set to zero "by hand" . While all terrestrial experiments are consistent with 
this assumption, it is now clear that the solution of the solar neutrino problem 
requires nonzero neutrino masses and neutrino mixing [1]. Another motivation 
for neutrino masses comes from cosmology with the need of a hot dark matter 
particle in the few eV mass range or of a heavy (~ MeV) unstable particle 
[2]. In both cases, the r- neutrino is an ideal candidate. However, astrophysics 
and cosmology provide also stringent limits for neutrino masses and lifetimes 
[3]. In particular, the requirement that the neutrino energy density does not 
overdose the Universe restricts the sum of the masses of all stable neutrino 
species to be less than 92Q u h 2 eV, where Q u h 2 < 1. 
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In extended SMs with neutrino mixing, neutrinos can decay by the process 
vu — > vl + 7 and, if the mass itlh of the heavy neutrino is higher than 
2m e , additionally by u H — > z/ L + e~ + e + . The resulting lifetimes r = 
are, although strongly model-dependent, extremely long. In the simplest case 
where right-handed singlet fields N m are added to the SM the decay rate of 
the resulting Dirac neutrinos is [3] 
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where m,H-,rnL are the masses of the heavy and of the light neutrino, respec- 
tively. If we assume that the r-lepton dominates the sum over the internal 
leptons / and that the elements of the mixing matrix U are of order unity, eq. 
(1) yields a lifetime 0(r SM ) = 10 29 years for m H = 30 eV and m L <C m H . 

Since tsm ^ ^Universe, neutrinos in the SM with lepton mixing can be treated 
as stable and the limit m v < 92 eV applies. The principal reason for the 
extremely small decay widths of neutrinos is their small mass. Obviously, the 
decay rate T for the radiative decay has to be proportional to aG 2 F , and since 
for rriH ^> the only energy scale available is the mass of the heavy neutrino, 
it follows T oc aGpm 5 H on purely dimensional grounds. 

This changes drastically if an external field is present. Then, in the low-field 
limit B <^ B CTti = mf/e the decay rates in magnetic fields B depend only on 
the dynamical field parameters 

Xi = m-^{p,eF^f = P —-^~, (2) 

where pj_ is the momentum perpendicular to B of the initial particle and B cr i 
are the critical fields of the charged particles, while in the strong field limit 
B ^> _B cr j = mf/e the decay rates can depend separately on p± and B [5]. In 
both cases, the essential point is that the energy scale is set by the charged 
particles and not by the mass of the decaying particle. 

The best example for this mechanism is photon splitting 7 — > 27. This process 
is forbidden in vacuum not only due to the Furry theorem but also due to the 
zero mass of the photon. In an external field photon splitting is allowed and 
the energy scale is set by the charged particles running in the loop. For, e.g., 
B m\je and uj 2m e the dominating contribution to the loop is that of 
the electron, and the decay rate is given by [4] 

w!^y7££) 8 . a ^. (3 ) 
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Consequently, the decay of all light particles should be enhanced in strong 
external fields. For the example of the radiative decay of neutrinos this was 
shown for various external fields configurations by Gvozdev et al. [6,7]. In 
the case of a weak magnetic field B <C B cr , they derived approximate decay 
rates valid for arbitrary momentum p± of the initial neutrino. However, in the 
strong field limit B B cr , they gave an asymptotic expression for the decay 
width valid only for p± <C 2m e [6]. 

Recently several authors proposed scenarios in which a strong magnetic field 
with B ^> m 2 e /e in the early Universe was created. Among the ideas considered 
are the production of a magnetic field of O(B) = 10~ s MQ VT /e « 10 43 G by the 
ferromagnetic Yang-Mills vacuum [8], of O(B) = rn^/e ~ 10 24 G during the 
electroweak phase transition [9], of O(B) = (A'Mpim0o) 2//3 ~ 10 48 G during 
hybrid inflation, where is the inflaton field with mass m and Higgs-coupling 
A' [10], inflation in string cosmology [11], etc. Subsequently, several authors 
used nucleosynthesis to derive bounds on primordial magnetic fields [12-14]. 
In these works, the effect of the magnetic field on the weak interaction rates 
n ^ p + e~ + v e and on the expansion rate of the Universe was taken into 
account but the neutrinos were treated as stable. The limits obtained vary 
between B « 2 x 10 13 G in ref. [13] and B > 1 x 10 15 G in ref. [14] at the 
beginning of nucleosynthesis (T « 1 MeV). 

Obviously, if the neutrino lifetime becomes of the order of the age of the 
Universe at the time of nucleosynthesis or before, they can no longer be treated 
as stable. Since the effect of a magnetic field is equivalent to an increase of the 
effective number of neutrino species N v , while the decay of a heavy neutrino 
before nucleosynthesis decreases N v from three to two, the limits obtained in 
ref. [13,14] could be even weakened by taking into account neutrino decays. 



2 Decay rate of v H — > v L + 7 



We consider the radiative decay of a heavy neutrino vh into a lighter neutrino 
vl an d a photon 7 in a strong magnetic field B 3> B cr = mj / e. Without loss 
of generality, we can choose the magnetic field as B = Be z and the three- 
momentum p of the heavy neutrino perpendicular to B, e.g. p^ = (E,p x , 0, 0). 
The momenta of the light neutrino and of the photon are denoted by p ,fl = 
(E',p^) and = (w,k); the components perpendicular to B have the index 
_L. If the energy of the heavy neutrino is small compared to the mass of the 
W-bosoia, E mw, the four-fermion interaction can be used and the decay 
is described by the Feynman diagram shown in Fig. 1. In this limit, Gvozdev 
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et al. [6] found for the matrix element Ai of the process vh — » vl + 7 

^(fi»fi cr ) = --^^C £ U Hl U* Ll e -^f(x), (4) 

V * l=e,[i,r I 

where an effective matrix element 

C = k\ (e^j v ) - (k^ Ju ) (k^el) (5) 

and a function f(x) depending only on the ratio x = 4mf jh\ were introduced. 
The other abbreviations have the following meaning: The neutrino current j 
is given by j M = z/ L (p 2 )7 M (l + ^)u H {p l ), = F^F^/B 2 - F^/B and F is 
the dual field tensor. For x > 1 the function / is given by 

3 

f(x) = -x (xfiarctanfi — 1) (6) 



with 



1 = = , (7) 
'4m 2 - k\ 



where we corrected a typographical error (in the argument of arctan) of ref. 
[6]; a check of the relative signs in eq. (6) is given by the requirement M. — > 
for uj — > 0. The function / can be analytically continued in the complex k\ 
plane. Then, for x < 1, the real part of / follows simply as 

3 

»(/) = -x (-x<£> arctanh $ - 1) (8) 

with 

$ = . ^ , (9) 
ykj_ — Am 2 

while the imaginary part of / is given by the dicontinuity across the branching 
cut starting at k\ = Am 2 , viz. 2i$s(f) = f(k\ + ie) — f(k\ — ie), 

3(/) = ^vrx 2 $. (10) 

Using arccothz = | In + ^ni for $$(z) < 0, one sees that f(x) for x < 1 
given in ref. [6] coincides — except for the sign of the term 3/2x — with the 
real part given by usQ but that the authors of ref. [6] missed to evaluate the 
imaginary part of /. 

In the following, we restrict ourselves to the case that the internal loop is 
dominated by the electron, i.e. that m 2 /e <C B <C rtv }Je. In this case, the 

2 This sign error has also been corrected by Gvozdev et al. in ref. [15]. 
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contribution of additional charged particles, e.g. from an enlarged Higgs sector, 
to the SM can be neglected. Therefore the decay rates will hold in the limit 
m,H 3> rriL n °t only for the SM with lepton mixing but also for more general 
models of massive neutrinos^} Moreover, we set UhiU^ = 1 for simplicity. 

The neutrino decay rate T is given by 

T = — *— r [ E dcu\M\ 2 . (11) 
Averaging \C\ 2 over spins, we obtain 



\C\' 2 = 16 kl(2Euj 2 E' - 2Euok z E' + Ek 2 z p' z + Ep' z u 2 



+2(e'uj - E'k z + p ' z uj + p' z k z ) (k z - iS)Euk z + Ep' z k] 
+ (EE'-p x p' x ) 
- k 2 x u 2 - k\k\ 

-2Eu 2 k 2 z E' + EuklE' + Ecu 3 k z p' 2 - Euk 3 z p' z 



k 2 y u 2 - k 2 y k 2 z + Eu 3 k z E' 



(12) 



In the limit E S> m^, m^, we can set E ~ p and E' ~ p' . Therefore the three 
particles propagate collinearly and the above expression reduces to 



\C\ 2 = 32u 4 EE' . 



(13) 



The differential decay rate dT/ du can be computed only for k\ ^ 2m e . 
Moreover, as the singularities of \f(x)\ 2 at k± = 2m e are not integrable, the 
total decay rate V above the pair creation threshold is ill-defined. Physically 
speaking, we have not taken into account the finite lifetime of neutrinos and 
photons in magnetic fields. Usually, the finite lifetime is incorporated by the 
replacement of the energies E by complexified energies E - \iT(E) in the 
denominators, where T(E) is the decay width of the state with energy E. We 
note that, although this replacement seems natural, it is not unambiguous. 
To proceed, we assume that the decay widths of the neutrinos are neglectable 
compared to those of photons and insert the total decay width 

r tot = r pair (^) + r splitt H (14) 

in the energy denominators of Q and $. Here, r pa i r denotes the decay width 
of the process 7 — ► 2e [16] and T sp \i tt the decay width of the process 7 — » 27 
[17]. 

Let us now consider the two analytically tractable limits p± <C 2m e and p± ^> 
3 If neutrinos are Majorana particles, the decay rates should be divided by two. 
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2m e . The function f(x) behaves as f(x) ~ | a; for a; — > and goes to 1 for 
rr — > oo. Inserting this and eq. (13) into eq. (11), we obtain 

r =i5^i( £sin ^ 6 (|;) 2 for ™««^«2">« (15) 

olG 

T = ^f^Xe for 2m e «p ± «m H ,, (16) 

— * 

where denotes the angle between the magnetic field B and p. The first ap- 
proximation was already derived in Ref. [6]. Note that the energy dependence 
of the total decay rate is the same as that for photon splitting. 

We now discuss our numerical results, which were all computed using the 
approximation eq. (13) and with 5 = 10 B CT , 9 = 90°. In Fig. 2 the differential 
decay rate dT/ doo (normalized to unity when integrated over the abscissa) 
is shown for energies of the initial neutrino below (left) and above (right) 
the pair-creation threshold as a function of the energy uj/E\ of the emitted 
photon. The energy distribution between the photon and the light neutrino 
becomes more and more asymmetric for E — ► 2m e , and finally, for E > 2m e , 
the resonance for uj = 2m e appears. For u = 100m e , the differential decay rate 
has a plateau similar to the case of photon splitting. 

A comparison between the total decay rate eq. (11), its low-energy approxi- 
mation eq. (15) and its high-energy approximation eq. (16) is made in Fig. 3. 
It can be seen that the exact rate is well reproduced by the approximations 
except for small region around 2m e . 

Finally, we have computed the lifetime r of the heavy neutrino as a function 
of the age t of the Universe. Hereby we have assumed that the magnetic field 
scales like B(t) = B(to)(to/t) 2 , and used as the starting value Bq the limits 
obtained in ref. [13] and [14] for the onset of nucleosynthesis at to ~ Is. If 
the lifetime is below t — r and the temperature T of the Universe lower than 
ran, the heavy neutrino decouples. Decoupling occurs for m H > 11 MeV using 
B = 2x 10 13 G and for 4MeV using B = 1 x 10 15 G. Therefore, the frequently 
discussed r-neutrino with m# ~ MeV [18] is possible in this scenario without 
introduction of new interactions. 



3 Summary 

We have presented the first analysis of the rate of radiative neutrino decay 
not restricted to momentum of the initial neutrino p± 2m e in the limit of 
strong magnetic fields B ^> B CT . We have derived two simple approximations 
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for the total decay rate valid for p± <C 2m e and p± 3> 2m e in this limit. The 
decay rates obtained are extremely enhanced by the strong magnetic field. In 
particular, if a strong primordial magnetic field existed in the early Universe, 
it is possible that a heavy neutrino with only standard model interactions may 
have decayed before the onset of nucleosynthesis. 
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Fig. 2. Differential decay rate dF/ du> (nor- 
malized to unity when integrated over the ab- 
scissa) as a function of u)/E\ for B = Wm^fe 
for (top) Ei/m e = 0.1 (full), 1.5 (dashed) and 
1.99 (dotted line), and (bottom) E\jm e = 2.1 
(full), 10 (dashed) and 100 (dotted line). 
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Fig. 1. Total decay rate T/s -1 for B = lOm^/e as a function of E\jm e (diamonds: 
exact rates, dashed straight lines: low-energy and high-energy approximation, re- 
spectively). 
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